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Abstract 



the subthreshold production of kaons and antikaons in heavy-ion collisions 



The in-medium mass and energy of kaons and antikaons are studied within 

m , 

the Relativistic Mean Field approach and compared with predictions from 

>. 

OO , chiral models by taking care of kaon-nucleon scattering data. Implications for 

in 
o 
t> 
o . 

\q , are discussed. We find only small corrections due to in-medium effects on the 

Os 

mean- field level for the relevant production processes for kaons. The produc- 
tion of kaons is even less favourable at high density due to repulsive vector 
interactions. We conclude that one has to go beyond mean- field approaches 



and take fluctuations and secondary production processes into account to 



explain the recently measured enhancement of kaon production at subthresh- 
old energies. The situation is different for antikaons where in-medium effects 
strongly enhances their production rates. We also see strong in-medium mod- 
ifications of the annihilation processes of antikaons and A's which might be 
visible in flow measurements. At high density, we predict that the threshold 
energy for antikaon and A production and annihilation become equal leading 
to similar numbers of antikaons and A's in the dense zone of a relativistic 
heavy ion collision. 
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I. INTRODUCTION 



In-medium properties of hadrons have received considerable attention recently, both 
experimentally and theoretically by studying relativistic heavy-ion collisions. Charged kaons 
(K + ) seem to be a quite promising tool for probing the dense interior of the collision zone 
as their mean free path is long enough to escape without further interactions. Kaplan and 
Nelson proposed first that a kaon condensed phase may be formed in the dense matter 
created in heavy ion collisions PJ. Further studies within the Nambu-Jona-Lasinio model 
0], chiral perturbation theory O] and an one-boson exchange model M showed that the 
kaon (K + ) sees a repulsive potential in the medium and will not condense. 

On the other side, the antikaon (K~) feels a strong attraction which is confirmed by recent 
calculations taking into account the contribution coming from the A(1405) resonance just 
below threshold ]5]J§. It was then predicted by chiral perturbation theory that a antikaon 
condensed phase will form in the dense interior of a neutron star |7| consistent with scattering 
data [Q,|J and Kaonic atoms 0. This approach has been criticised in PJlO[| as the scalar 
density is set equal to the baryon density and higher order terms in density are neglected. 
The appearance of hyperons Jll| shifts also the onset of a condensed phase to higher density. 
As shown in [0 a strong nonlinear dependence on density and the implementation of 
hyperon-hyperon interactions even prevents an antikaon condensed phase inside a neutron 
star. 

Here we will continue our work for neutron stars [1^] and apply it for the situation in 
heavy ion collisions at threshold. Subthreshold production rates of K + in heavy-ion reactions 



were recently measured at GSI [p^| . Earlier work showed that the in-medium modifications of 
kaons and antikaons might be measurable in heavy ion collisions at threshold. For example, 
it was shown that the kaons are sensitive to the equation of state (EOS) ||14|-[T6[|. A softer 
EOS produces more kaons than a hard one. On the other hand it will also depend on the 



parametrisation used for the cross section [17,18] but not on N-body collisions p9l and not 



on the high-momentum tail of the nucleons [20|. The influence of rescattering and formation 
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of resonance (A) matter was studied in the QMD model and the RBUU model [21| and 
it was demonstrated that they are essential to explain the data. In-medium modifications 
of the effective energy of the kaon were studied in p2[ using again the RBUU model. The 



results are essential similar to the ones obtained without medium modifications ||21|| , because 
the in-medium kaon mass used is quite close to the respective vacuum mass. But there exist 
other observables which might be better suited for extracting in-medium effects. The flow of 



kaons might be a promising tool for measuring the kaon potential in dense matter ||23|| . And 



more pronounced in-medium effects are expected for the case of K~ |24]]. Indeed, enhanced 
production rates for K~ have been seen at GSI recently ||25|| . 

In this paper we want to examine the possible influence of a dense nuclear environment 
on the properties of kaons and antikaons. We show that the in-medium effects on the mean- 
field level can not explain the measured enhanced production rates of kaons in contrary to 



the conclusion drawn in ref. [p4}| . We discuss two different approaches: first an one-boson 
exchange model and second a chiral approach where the parameters are fixed by s-wave 
scattering lengths and the low density theorem. In- medium effects for A's are also taken 
into account by linking them to hypernuclear data. We show that the phase space in the 
medium does not change considerably for the processes NN— >NAK and secondary processes 
as 7rN— > AK and NA ^NAK due to cancelation effects. On the contrary, effects nonlinear in 
density even cause an enhanced repulsion at highly dense matter for these processes. Hence, 
subthreshold production of kaons seems not to probe the potentials of the very dense region 
of a heavy-ion collision. On the other hand, in-medium effects are essential for explaining 
the enhanced production of antikaons. We show that the process NN— >NNKK is enhanced 
in the medium while the annihilation process XN— » Air is reduced. We also find that the 
annihilation process AN— ►NNA' is essentially enhanced in dense matter and might be equally 
important as antikaon annihilation. This behaviour may lead to equal numbers of antikaons 
and A's in the dense zone of a relativistic heavy ion collision. 

The paper is organized as follows: first we introduce the Relativistic Mean Field (RMF) 
model and extend it to include A's. In the second section we discuss two different interaction 



schemes for the kaons with nuclear matter, one based on an one-boson exchange model and 
the other on chiral perturbation theory (ChPT). The parameters are fixed to the s-wave KN 
scattering lengths. Results for the in-medium effects on kaon and antikaon production are 
presented in the third section. The last section is devoted to conclusions and an outlook. 



II. THE RMF MODEL 

The RMF model has been proven to give a good description of nuclear matter in bulk 
and of the properties of finite nuclei [ 2"6 , |27J . We start from the Lagrangian 



£ = ^(iyfy - m N )* N + -d^ad^a - U(a) 



-\g^G, v + ifn* V% - -B^B, V + \m 2 p R?R, 



(1) 



where the nucleons interact via an attractive scalar (a) and a repulsive vector (V^) meson 
field. The term U (er) stands for the scalar selfinteraction 



U(a) 



1 9 9 b Q C A 

-mia + -cr + -a 

2 a 3 4 



(2) 



introduced by Boguta and Bodmer [28 to get a correct compressibility of nuclear matter 



(for another stabilized functional form see [29]). The parameters of this Lagrangian can be 
fixed to bulk properties [28] or to the properties of finite nuclei [30,j3l|. A general discussion 



about the scalar selfinteraction terms can be found in [32]. Bodmer proposed an additional 



selfinteraction term for the vector field [33] 



c V 4 = -d(v„v»y 



(3) 



4 -vA. 

which leads to a soft equation of state at high densities in agreement with Dirac-Briickner 
calculations P3,B3. Fits to the properties of nuclei with this new term are quite successful 



| 3q| . In the following we take mostly the parameter sets NL-Z [3(| which is the commonly 



used parameter set NL1 with a better zero-point energy correction and the recent set TM1 
p5| with vector selfinteraction terms. The former one gives a rather stiff equation of state 
while the latter one a rather soft one. 

The implementation of hyperons proceeds as 

C A = ^ h {iYd^ - m A )^ A - g aA ^ A y A a - ^a^aT^a^ (4) 

and the two new coupling constants can be fixed to hypernuclear data |36||. The main feature 



of hypernuclei is that the depth of the A-potential is about 

U ( A N) = g aA a^ + g uK V^' « -30 MeV (5) 
in saturated nuclear matter which already fixes one coupling constant of the A, say g aA 



37,38]. The vector coupling constant g ujA is then given by 



9wh = ^gu>N (6) 

when using SU(6)-symmetry (the quark model, see e.g. |^|). The SU(6)-symmetry also 
secures that the spin-orbit force is negligible small as there is no experimental evidence for 



a spin-orbit splitting for hypernuclear levels. Noble showed first [^0[ that the contribution 
of the vector terms to the spin-orbit term nearly cancel each other when taking into account 
the tensor force and SU(6)-symmetry. The tensor force vanishes in bulk matter on the mean 
field level as it is proportional to the gradient of the fields. Therefore it is not considered 
here. 

The in-medium energy of nucleons and hyperons is then given by 



E N {p) = J (m N + g aN a) 2 + p 2 + g^Vo + g P NT Ro,o (7) 



Ea(p) = V (m A + gainer) 2 + p 2 + guhVo ■ (8) 

It is important to note that the parameters here are connected to properties at normal 
nuclear matter density. The in-medium effects for nucleons and A's at this point are known 
and should be taken into account when studying the influences and the signals of a dense 



nuclear environment. As pointed out in j|T] three-body forces are also important to explain 



hypernuclear data. As these forces are repulsive, the hyperon potential shows a nonlinear 
behaviour with density and changes sign at higher density. Figure |1] shows the Schrodinger 
equivalent potential defined as 

^sev = g*A° + 9uaV + ((s^aoO 2 - (fiU\Vb) 2 ) (9) 
for different parameter sets in comparison with the findings of the non-relativistic approach 



4341 . The overall behaviour is quite similar despite of the EOS used. The nonlinear behaviour 
of the scalar field with density seems to simulate the repulsive three-body force of the 
nonrelativistic approach. It also demonstrate that is crucial to make a difference between 
scalar and vector (baryon) density. This turning of the hyperon potential will be quite 
important for our discussion of the kaon production in the medium. 

III. KAON INTERACTIONS 

The case for the kaon is quite distinct from that of the A. There does not exist any 
kaon-nuclear states similar to hypernuclei as the KN-interaction is known to be repulsive. 
Taking the (real) isospin averaged KN-scattering length qkn = (3ag =1 + ag =0 )/4 = —0.255 



fm [42| and using the low density theorem one gets a repulsive optical potential depth at 



normal nuclear density of about 



U™ t = -— (l + —) a KN p N « +29 MeV p N /p (10) 



compatible with kaon (K + ) scattering on nuclear targets [|39" |. Here we have taken the 
groundstate density to be p = 0.15 fm~ 3 . The repulsive interaction is the reason why 
kaons have a long mean-free path in nuclear matter. Note that the potential depth is 
just opposite to the one of the A which signals a significant cancellation of attractive and 
repulsive terms in the medium. On the other hand, a recent experiment measured an 
enhanced cross section for K + scattering on nuclear targets incompatible with multiple 
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scattering arguments. This is the so called K + -puzzle which is still unresolved. The isospin 
dependent potential in nuclear matter can be estimated from the isospin scattering length 
fliso = ( a o =1 — a o = °)/4 = -0.055 fm and the low density theorem 



rKN 27F ft , m K \ - 



= 1 + — «isoPiso « +6 MeV Piso /p (11) 

rriK V m N J 

where pi SO is the isospin density of the system. For led one can estimate pj SO ~ (2Z — 
A)/ApN ~ — 0.21pAr, which gives about 1 MeV correction at normal nuclear density. 
For antikaons the annihilation processes 

A > + N->Y + vr (Y = A,E) (12) 

gives a big imaginary part for the scattering lengths. At first glance the experimental 
situation seems to be contradictory: The available K~N-scattering indicates a repulsive 
interaction while the K~-atomic data demands an attractive potential. The situation can 
be remedied by taking care of the existence of the A(1405)-resonance just below threshold. 
Recently an improved fit of K~-atomic data was carried out assuming a nonlinear density 
dependence of the effective t- matrix Ell . It has been shown that the real part of the antikaon 



optical potential can be as attractive as 

« -200 ± 20 MeV (13) 

at normal nuclear matter density while being slightly repulsive at very low densities in ac- 
cordance with K~p-scattering. The change of the sign and the nonlinear density dependance 
results from the A(1405)-resonance. Also another family of solutions have been found with 
a moderate potential depth around —50 MeV. Note that also the standard linear extrapo- 
lation gives only values of about —85 MeV JPJ. These latter two solutions are not getting 
repulsive at low densities, i.e. fulfilling the low-density theorem. 

The K~N-scattering data can be explained by vector meson exchange models where the 
A(1405) is a quasi-bound state in the t-channel PHHBJ. In a recent paper the coupled channel 



analysis of Siegel and Weise |45j has been also applied for interactions terms coming from 



chiral perturbation theory The coupled channel formalism automatically generates the 
A(1405) and successfully describes the low energy K~p-scattering data. 

In the following we adopt the meson-exchange picture and the chiral approach for the 
KN-interaction on the mean-field level and fix the parameters to the KN-scattering length. 
The case for the antikaons is then given by a G-parity transformation which simply changes 
the sign of the vector potential term. This simple treatment does not take care of the 
important contribution of the A(1405) resonance. But there exist some hints that this 
resonance seems to be less important in dense matter (which happens when the antikaon 
energy is shifted down below m(A(1405)) — m^ « 466 MeV). In ref. || a separable potential 
was applied for the K~p-interaction for finite density. Indeed, it was found that the mass of 
the A(1405) is shifted upwards and exceeds the K~p threshold already at densities of about 
p ~ 0.4p - This is supported by recent findings within a chiral approach ||, where this 
resonances vanishes at very low densities p ~ 0.2po due to Pauli-blocking effects. In this 
case the use of mean-field potentials may be justified. Hence, we simplify our calculation by 
neglecting the contributions coming from the A(1405) in the medium and treat the problem 
on the tree-level using G-parity. Nevertheless, the results presented for the NN— »NNKK 
case should be taken with some care. More elaborated models are needed to draw final 
conclusions about the in-medium property of antikaons in the medium. 



A. One boson-exchange approach 

We start from the following Lagrangian [|12|] 

C KN = D^KD^K - m\KK - g aK m K KKo - g &K m K KrKb (14) 

with the covariant derivative 

D^ = d^ + ig^xV^ + ig^rR^ . (15) 

For completeness we also add isospin-dependent terms which couple to an isovector-scalar 
field (6) and an isovector- vector field (R/j). Note that interaction terms of the form 

8 



Ckna = -9kna (Nr l5 AK + Kt^NK) (16) 

do not contribute on the mean-field level as they are off-diagonal terms. We will come to 
this point later in more detail. 

The coupling constants to the vector mesons are chosen from the SU(3)-relations assum- 
ing ideal mixing 

2sw = 2g pK = grrnp = 6.04 (17) 

where g nnp is fixed by the p decay width. The scalar coupling constants can be fixed to 
the s-wave KN-scattering lengths fl2f . The isospin averaged scattering length in the tree 



approximation is given by 



a KN = 1" + = . „ /"* ■ , ^ - 2™ = "0.255 fm (18) 

4 4 47r (1 + mK/TriNj \ m l a m£ J 

where only the isoscalar terms contribute. This can be used to fix g a x for known g w x = 3.02. 
The KN-scattering lengths for a given Isospin / on the tree level are then given by 



a, 



1=1 rn K 



( gaKgaN _^ gsKgsN _ ^gvKgujN_ _ ^gpKgpN_\ 

11 An (1 + ttik /ttin) \ m a m 1 m t m p / 

a I=0 _ IJ^K ( gaKgaN _ ^ SSKdSN _ ^9wK9ujN_ ^9pK9pN_\ ^qn 

Air (1 + m^/m^v) y m 2 m 2 m 2 p J 



Recent experimental values are a ! 1 = 0.31 fm and = —0.09 fm [42]. The importance 



of the 5-meson exchange contribution can be seen by looking at the a 1 ^ scattering length. 
The vector terms largely cancel each other as g^Kg^N ~ 3g p xg P N- Hence, without the 
contribution from the ^-exchange one gets 

af° = , mK , : ( 9 -^] « 0.4 fm (21) 
4tt (1 + m K /m N ) \ J V ' 

in contradiction with experiment (here we used g a ^ = 10, g^N — 13 as standard values for 
the RMF model). Including the £-meson term and using g^N = 5.95 from the Bonn model 
| 49| one can fit both scattering lengths nicely for 



g aK « 1.9 - 2.3 , g SK « 5.6 - 6.4 (22) 



for the various nucleonic parameter sets used in the literature (see Table |). Note that the 
values of g a x significantly deviate from the simple quark-model (simple quark counting gives 
9aK = 9aN /3 ~ 3.3). The coupling of the kaon to the 5-meson is quite strong. Therefore, 
we expect some effects for isospin-asymmetric systems which we will discuss later. 

The fit based on the adjustment to the KN-scattering lengths leads to an optical potential 

of 

2m K U* t = u% - m 2 K = g aK am K - 2g ujK uj R V - {g^Vo) 2 (23) 
which gives U^ t = —(85 -r 100) MeV at normal nuclear density for the parameter sets 



used. These values are lower than the ones quoted in our previous work [12|] as we use here 



the vacuum kaon mass nix instead of the reduced mass hkn- We think that this is more 



consistent with the parametrisation used in the study of Kaonic atoms [44], but now our 



value is much closer to the standard fit which gives U* t = -85 MeV. Note that the optical 
potential as defined in ( p3| ) is always lower than the relativistic potential the kaon feels at 
normal nuclear density which is about 

U* h =u R -m K ^ -(95 110) MeV . (24) 



This definition corresponds to the sum of scalar and vector potentials as discussed in [pO 
Nevertheless, the scalar and also the vector potential are much lower than the ones deduced 
from simple quark model counting as used in |5(J . The reason is that our coupling ratios are 
about g a K/g<rN ~ 1/5 and g^x/g^N ~ 0.23 (see Table |) which significantly deviates from 
the simple quark model value of 1/3. 

We have also studied the influences of off-shell terms which have only small influences 
on the in- medium behaviour of kaons (see [0). Note that off-shell terms are not needed for 
describing the s-wave KN-scattering lengths correctly. On the other hand, they are essential 
for the chiral approach which we will discuss in the next section. 

The equation of motion for kaons in the mean-field approximation in uniform matter 
reads 

[d^ + m\ + g aK m K a + gsKm K T 5 + 2(g wK V + g P KToRo fi )id^ 
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-(9*kV + g pK T Ro,o) 2 } K = . (25) 

Note that there appears terms quadratic in the vector fields in eq. (p5|). The importance of 
the isospin dependent terms can be estimated from the equation of motions for the vector 
fields in uniform matter 

mlV + dV* = g„ N p N (26) 
m 2 p Ro >0 = 9 p n{p p - Pn) (27) 

where p p and p n are the densities of protons and neutrons, respectively. For led one gets 
P P — Pn ~ (2Z — A)/ApN » — 0.21pAr, and the isovector correction is then about 

« 21%^| « 8% 

if one neglects the vector field selfinteraction which holds for low densities. Hence, the 
isovector contributions are expected to be small for the densities considered here (p < 3po)- 
The effective mass of the kaon is given by 



m *K = ym 2 K + m K {g aK a + g SK r S) . (28) 

The scalar field a reduces the effective mass of the kaon in the medium, i.e. the scalar 
interaction is attractive. The isovector-scalar field 5 shifts the effective mass if there is 
an isospin asymmetry in the system. Note that for kaons as bosons the dependence on the 
scalar potential is different from that for nucleons (fermions) : for low densities the reduction 
of the kaon mass in the medium is proportional to the square root of the scalar attraction 
while it is linear for the case of baryons (see eq. (|7P). Moreover we point out that the scalar 
potentials always follow the scalar density as demanded by Lorentz invariance. The scalar 
density is saturating in dense matter to ensure the existence of a saturation point of the 
equation of state. As shown in [|],|n| these nonlinear effects are important already at a 



moderate density and causes a saturation of the effective kaon mass with density. We can 
even go further and say that there exists a minimum effective kaon mass. As th effective 
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mass of the nucleoli approaches zero at high density {m* N — > 0) in the Walecka model, the 
minimum scalar field is about a min = —rri]si/g C N and one gets for the minimum effective 
kaon mass 



m* Kmin = Jm 2 K + m K g aK a min « Jm 2 K + mK9aK ~ 390 MeV (29) 

V V m N9aN 

for the parameters of Table |. 

Decomposing the kaon field into plane waves one obtains the following dispersion relation 
for kaons (upper sign) and antikaons (lower sign) in uniform matter composed of nucleons 
only 



^k,k = y m K 2 + k 2 ± (guxVo + g P KT Ro,o) ■ (30) 

Note that due to the covariant derivative coupling scheme flT5|) the vector term appears 
linearly in the kaon energy. The vector field is repulsive (attractive) for the kaon (antikaon) 
and will dominate the behaviour in very dense matter. For high density the kaon (antikaon) 
energy is then increasing (decreasing) as p 1//3 because the vector field is growing with p 1 / 3 



if one takes into account the vector field selfinteraction term (see eq. (26)). Otherwise it is 
changing linear in density. 

B. Chiral Approach 

We follow the procedure outlined in || starting from the Nonlinear Chiral Lagrangian 
in next-to-leading order 

3i 



^chiral 

KN — ~~gp 



K 



N'j^N (Kd^K) + Nr-f^N (Krd^K 



v C - 

~^ — j^-NNKK + -^NtNKtK 
Ik Ik 

+^NN (d^Kd^K) + ^NtN (d^Kfd^K) (31) 
Ik Ik 

where fx = 93 MeV is the kaon decay constant and S^tv is the KN sigma term. The first 
two terms are the Tomozawa- Weinberg terms and are in leading order of the chiral expan- 
sion. These are vector interactions terms and repulsive (attractive) for kaons (antikaons). 
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The other terms are in next-to-leading order. The next two terms are scalar interactions 
which will decrease the effective mass of the kaon and antikaon. The last two terms are 
the so called off-shell terms which will modify the scalar attraction. Here one encounters 
striking similarity with the RMF model as the interaction is governed by scalar and vec- 
tor interactions (see [50| for a discussion about this point). In the original paper || the 



authors choose T*kn ~ 2m n in accordance with the Bonn model |]46 |. More recently the 



value T*kn = 450 ± 30 MeV is favoured according to lattice gauge calculations |5l| . The 



constant C can be fixed from the Gell-Mann-Okubo mass formula to C = 33.5 MeV. The 
corresponding scattering lengths 

afi? = A n n r 7 V \~ m K + ^kn + C + (£> + D') m\ 

af-o = A J2 m 1 / V \+Zkn -3C+(d~ 3D') m 2 K ] (32) 

4tt/& (1 + m K /m N ) L V J K \ 

determine the constants D and D' for a given T<kn via the relations 

D « 0.33/m^ - T, KN /m 2 K , & « 0.16/m* - C/m 2 K . (33) 

Note that the off-shell terms involving the constants D and D' are essential for a correct 
description of the scattering lengths (see || for details). The equation of motion in the 
mean-field approximation and in uniform matter reads 



I HO 



( a an i 2 ^KN C 

+^-p s ^cP + W T ^s°d^ + ^PNdt + -^r oP %°d t ) K = (34) 
jk Jk ^jk ^jk J 

where p l s so = p SjP — p S) „ is the scalar-isovector density and p l ^° = p p — p n the vector-isovector 
density which are simply the difference of the corresponding densities of protons and neu- 
trons. The mass of the kaon is shifted by 



m K = \ m K - -J^Ps ~ -TTToff (35) 

V JK JK 



and the same arguments as for the case of the one-boson exchange model holds. One sees 
again that the scalar potential for the kaon behaves differently as the one for nucleons (0). 
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More important is, that there exists a minimum effective kaon mass as a minimum scalar 
field implies a maximum scalar density for the Walecka model which is about p s , m ax ~ 2/?o 
(see e.g. |[29|| ). This gives a minimum kaon effective mass of 350 — 400 MeV depending on 
the kaon-nucleon sigma term. The influence of the isovector terms can be estimated from 
isospin considerations: e.g. for led one has p p — p n ~ (2Z — A)/Ap N w —0.21p Nl i-e. about 
21/3 = 7% correction for the vector-isovector term of eq. (|3~4"1) . This is in accordance with 
our estimate for the one-boson exchange model in the previous section. In the following we 
will neglect the isovector contributions. Fourier transformation of the equation of motion 
yields 

- u 2 + k 2 + U(u, k; p N ) = -uj 2 + k 2 + m 2 K - —^p s - ^Ps^ 2 - tt^Pv = . (36) 

JK JK 4jK 

where H(uj, k; p^) is kaon self energy which depends in general also on the kaon energy. This 
has to be taken into account to get the energy of a kaon/antikaon in the nuclear medium 

-i 



^K,K 



1 + Jit>A • < 37 ) 



where m* K is defined by eq. (^1). Here we note that in the high density limit the kaon 
energy is growing linear with density while for the antikaon the energy saturates at m* K as 
the vector contributions cancel each other. The optical potential at normal nuclear density 

U*t = 7^—n(uj Rj k = 0; p ) « -68 MeV (38) 

is rather moderate while the relativistic potential is about —75 MeV. This is in contrast to 
the findings of Brown and Rho |5(J who gets a rather deep potential of —200 MeV. There are 
several reasons for this discrepancy: first BR-scaling is not taken into account here (which 
gives an additional factor of 5/3 at po, i.e. a potential of —125 MeV), second Brown and Rho 
neglect the off-shell term and do not take into account the KN-scattering lengths, third they 
neglect the energy dependence of the kaon self energy, fourth they assume that the scalar 
and vector density are equal, fifth they neglect that the scalar potential of the kaon behaves 
differently in matter compared to the nucleon one (see the discussion of eq. (|35|) above). 
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IV. RESULTS 



A. Kaon energy in matter 

Recently, the dynamics of the A(1405) has been studied in nuclear matter using a coupled 
channel formalism |5|||]. The most important finding is that the effects coming from the 
A(1405) vanishes at rather low densities (p < 0.25po)- The optical potential for the antikaon 
is about —100 MeV || and —107 MeV || corresponding to a kaon energy of uj^ = 380 MeV 
and lok = 372 MeV at normal nuclear density, respectively. These values are in accordance 
with the ones calculated in the mean field approximation in the previous sections within the 
relativistic mean field (RMF) model and the chiral perturbation theory (ChPT). 

In the following we discuss the in-medium energy of kaons and antikaons in nuclear 
matter using a soft (set TM1 of Table |) and a hard (set NL-Z of Table H) equation of state. 
Figure H shows the energy of kaons (upper curves) and antikaons (lower curves) with the soft 
EOS for the RMF model (eq. (0)), ChPT (eq. (p7|)) and the results of the coupled channel 
analysis of Waas et al. || . In the case of ChPT we discuss three cases: i) for a sigma term 
of Yikn = 2m„- as used in |J, ii) for a sigma term of Y*kn = 450 MeV as derived from recent 



lattice data [51], iii) for vanishing off-shell terms (denoted as D = 0) and a sigma term of 
Exjv = 2m n as used as input for the RBUU calculations p!6| , p2] -p4| . 

All models show a quite similar behaviour in Fig. [I] for the kaon energy at low density 
except for the case D — 0. This results from the low density theorem and is a generic feature 
when the coupling constants are fixed to the KN-scattering lengths. Neglecting the off-shell 
terms, i.e. setting D — 0, violates the low density theorem. This gives a slower raise of 
the kaon energy with density and the kaon energy nearly stays constant for a wide range of 
density. Note that this latter parametrisation for the kaon energy is used in the dynamical 



calculations p2[ . For higher density the other curves also start to deviate. The ChPT gives 
a higher kaon mass, i.e. more repulsion than in the RMF model. The results of the coupled 
channel calculation seems to follow more closely the one of ChPT. At p = 3po the kaon 



15 



energy reads ojk = 585 MeV for the RMF model, ojk = 630 MeV for the coupled channel 
analysis and u K = 640 — 670 MeV for ChPT, so they deviate about 85 MeV from each 
other. 

The antikaon energy (lower curves) of the different models is always attractive, except for 
the small density region for the coupled channel calculations due to the A(1405) resonance. 
The latter one gives the most attraction of about = 217 MeV at p = 3po, followed by 
the RMF model with uj r = 263 MeV and the ChPT with uj r = 280 4 300 MeV. The curve 
for the case of D = used in |24| follows closely the one for the RMF model. All the curves 
for the antikaon energy are lying surprisingly close together. Note that the prediction of 
ChPT is rather insensitive to the choice of T^kn but rather sensitive to the off-shell terms, 
especially for the kaon energy. 

In Fig. [3] the case for the hard EOS is plotted. Now the curves of the kaon energy are 
lying very close together, even at higher densities. This is due to the fact that the vector 
potential in the RMF model is now raising linear with density as in the ChPT in contrary 
to the soft EOS where it raises like p 1 ^ 3 due to the vector self-interaction terms. The energy 
of the kaon is now between uj^ = 630 4- 670 MeV at p = 3po- Without the off-shell terms, 
the kaon energy significantly deviates from the other curves and stays rather constant up to 
1.5po- Note that the overall changes for the hard EOS compared to the soft EOS are quite 
moderate, especially when using ChPT, and only show up at higher density. 

The different predictions for the antikaon energy seems to split now into two regimes: 
the results for the ChPT give a antikaon energy of about Uk = 300 MeV at p = 3po rather 
independent of the off-shell term and the choice of T^kn, while the RMF model and the 
coupled channel analysis get around = 200 MeV at p = 3po- The antikaon energy within 
the RMF model is now much deeper due to the stronger vector potential compared to the 
soft EOS. We want to point out again, that Dirac-Briickner calculations seems to favour the 



soft EOS p3| , |34|| . Nevertheless, we see that the differences in the kaon/antikaon energy due 
to the EOS are well within the differences of the model predictions. 
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B. Threshold energy for kaon production in matter 



In the following we discuss the shift of the threshold energy of various processes for heavy 
ion collisions due to medium modifications. 

Kaons are mainly produced at threshold via the process NN— >NAK. The minimum energy 
needed is Q(NAK) « 671 MeV in vacuum. In the medium, the threshold is shifted to 

Q(NAK) = E A (p = 0) + oo K (k = 0) - E N (p = 0) (39) 

where we assume that the outgoing nucleon is not Pauli-blocked in the hot zone of the 
collision. Hence, the subthreshold production of kaons is sensitive to three different in- 
medium effects: the EOS (En), the A potential (E\) and the kaon energy (ojk) in medium. 
These effects will partly cancel each other as the kaon feels a repulsive potential of 29 MeV 
(eq. (|10D ) while the A sees an attractive potential of —30 MeV at p (eq. (§)). Therefore, 
subthreshold kaon production seems to probe mainly the EOS. As the nucleons feel an 
attractive potential of about —60 MeV the threshold will be shifted upwards at normal 
nuclear density by this amount and the production of kaons is reduced in the medium. This 
is indeed the case as can be seen from Fig. f| which shows the threshold energy Q(NAK) as a 
function of density. The similar behaviour of the different curves at low density is due to the 
low-density theorem. At p = 3po the value of Q(NAK) reaches about 800 MeV for the RMF 
model and about 860 MeV for ChPT which is quite insensitive to the value of the sigma 
term. Without off-shell terms, the threshold energy is underestimated in medium, and we 



expect that the production rates for kaons calculated in P2| , ^4| are overestimated. Note, that 
all calculations ignoring in-medium effects |15|J21| will also give a too high production rate 
for kaons. The case for the hard EOS is plotted in Fig. |^. The behaviour of the threshold 
energy in medium is quite similar for the different EOS considered here. Again, the low 
density limit more or less fixes the shape of the curves of the kaon energy independent of the 
EOS used. The curves for the RMF model and ChPT are lying closely between 800 — 830 
MeV at p = 3p - Especially, the curve for the RMF model does not change considerable 
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for the hard EOS compared to the soft one. All curves seem to saturate for the hard EOS 
but are lying within the uncertainties of the different models used for the kaon energy. A 
definite conclusion whether or not subthreshold kaon production probes the EOS can not 
be drawn until the in-medium properties of the kaon can be determined more precisely. 

Antikaons are created in heavy ion collisions first by the process NN— s-NNKA. The 
threshold value of Q(KK) w 988 MeV is modified in the medium by the sum of the kaon 
and antikaon energy Q(KK) = ux(k = 0) + uj R {k = 0). Therefore, subthreshold antikaon 
production probes the in-medium property of kaons and antikaons solely. As the vector 
potential cancels out approximately, it will mainly depend on the scalar potential the kaon 
feels in the medium. The upper curves in Fig. f| show that indeed Q(KK) is reduced in 
the medium in all models discussed here. ChPT predicts an in-medium reduction of about 
—56 MeV at maximum compared to the vacuum and then the curves go up again for higher 
density. The reason is that the sum of the kaon and antikaon energy contains a term coming 
form the Tomozawa- Weinberg term 



u K (k = 0) + u R (k = 0) 



m *K (l + ^Ps) + ( 7^2 Pn) (l + ^rPs) , (40) 



which is repulsive and dominates at higher density. On the other side, the RMF model gives 
a reduction of about —140 MeV at p = 3po- The Q- value is steadily decreasing as the sum 
of the kaon and antikaon energy 



u K (k = 0) + uj R (k = 0) = 2m* K = 2^m 2 K + m K g aK a (41) 

depends on the attractive scalar potential only. The curve used in RBUU calculations with 
a soft EOS [f24| is lying even lower and hence, the production rates of antikaons seems to 
be overestimated. Using the hard EOS (Fig. [5]) the situation does not change significantly. 
The Q-value in the RMF model is now reduced by —160 MeV at p = 3po- The curves for 
the ChPT go up stronger at higher density compared to the soft EOS as they are sensitive 
to the strength of the vector potential (i.e. to the behaviour of the EOS at high density) in 
contrast to the RMF model. 



As an interesting fact, the Q- values for kaon and antikaon production are lying close 
together for the RMF model. Note that this does not mean that the numbers of produced 
kaons and antikaons are the same inside the dense medium. The kaons will be produced 
at different density and the average Q- value over the density profile will give a measure for 
the produced kaons and antikaons in the medium. On the other hand, the production of 
kaons will be dominated by the secondary processes (rescattering effects) NA — >NAK and 
7rN— ► AK, the production of antikaons by the processes NA — »NNK.K" and 7rN— >KK. Let us 
assume that the change of the A mass and energy is equal to that of the nucleon. Then the 
Q-values of these channels can be simply derived by shifting the corresponding curves for 
the Q- values of Figs. [| and |5] down by — ~ —290 MeV (ignoring the finite width of 
the A) and by —m^, respectively. If the A feels a higher (lower) potential than the nucleon, 
then this will suppress (enhance) subthreshold kaon production. Processes involving two 
A's in the entrance channel will decrease the Q-value by —580 MeV compared to the two 
nucleon one and hence, enhanced production of kaons will be sensitive to A matter (density 
isomers) [JTJJ. 

Also annihilation processes will play a dominant role at high density. Kaons will not 
annihilate and escape due to their long free mean path. But the charge exchange reaction 
K + n — > K°p will act like an annihilation process for kaons as only charged particles are 



measured in the present heavy ion experiments |]T3] . This process will be modified in the 



medium only by isovector potentials. We do not expect changes of the threshold energy for 



isospin symmetric systems. As the isospin potential for the kaon is negligible (see eq. [II]) 
the threshold energy will be only shifted by the isovector potential of the nucleons. The 
maximum effect will be seen for systems like led where one gets 

2 

Q(iso) « E n (p = 0) - E p (p = 0) = -2£ piVJ R 0i o = -2^f p iso » (16 - 19) MeV^ (42) 

m p Po 

with the parameters of Table |. Hence, the charge exchange process will be a little bit 
suppressed in isospin asymmetric systems. The change is quite moderate but comparable 
with the in-medium shift of the Q-value for the kaon production process. 
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Antikaons will annihilate strongly due to the process ^N— > An which is exothermal in 
vacuum (Q(KN— > An) w —180 MeV). Also the A can annihilate via the process An— >NNA", 
but this process is endothermal (Q(AN— >NNA^) ?s 317 MeV) in vacuum and is usually 
neglected. Nevertheless, we expect rather strong in-medium modifications of these Q-values 
as the antikaon energy is involved which changes considerably in nuclear matter. The lower 
curves in Fig. f| and |5| show these Q-values as a function of density. All models give an 
astonishingly similar strong behaviour in dense matter: Q(KN—> An) is going up with 
density and crosses zero at p ~ 1.5po while Q(AN— >NNK) is decreasing rapidly. For the soft 
EOS, both Q- values reaches even similar values at high density of about 70 MeV regardless 
of the model used. This means that the annihilation of A's is favoured in the medium while 
the annihilation of antikaons is suppressed. At very high density these processes are even 
equally possible. For the hard EOS (Fig. ^) the Q-values for the ChPT seem to saturate at 
high density at Q(KN^ An) « MeV and Q(AN^NN^) « 140 MeV. On the contrary, 
the curves for the RMF model show a crossing, so that the situation is reversed and one gets 
Q(KN^ An) « 115 MeV and Q(AN^NNA > ) rj 20 MeV. It would be interesting to examine 
how these in-medium effects of the annihilation process will influence the antikaon and A 
spectra in heavy ion collisions at subthreshold energy where they are most pronounced. 



V. CONCLUSIONS AND OUTLOOK 

The study of the in-medium properties of kaons and antikaons shows that it is important 
to link the models to the available data, here to the KN scattering data, and to take into 
account effects nonlinear in density. Then one gets rather similar predictions in the models 
discussed here for the energy of kaons and antikaons in nuclear matter up to a certain 
density, say (1 — 2)po- The kaon energy in matter is well determined by the low density 
theorem, while the antikaon energy is more model dependent. The threshold energy for the 
production of kaons is shifted up in dense matter, while the one for antikaons is considerable 
decreased. Also the threshold energy for the annihilation processes for antikaons and A's 
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show strong in-medium modification and can even get similar values at high density. Hence, 
it will be important to study the process AN— >NNK" in the medium which will enhance the 
number of produced antikaons in heavy ion collisions at subthreshold energy. This will also 
change the flow pattern of antikaons and A's and will cause e.g. an antiflow of A's for central 
rapidities. As the threshold energy for antikaon and A production as well as for annihilation 
become equal around p ~ 3po the number of antikaons and A's will be predicted to be 
equal in the dense zone of a heavy ion collision at subthreshold energy due to in-medium 
effects. The number of kaons will then be twice the number of antikaons due to strangeness 
conservation. This effect might be seen at midrapidity and high momenta. 

Insofar, we have only discussed effects on the mean-field level which will cause shifts of 
the threshold energy and essentially modify the phase space of the reactions in the medium. 
Using Fermis golden rule, one can now implement these modifications into a dynamical model 
by simply changing the energy of all hadrons consistently and leaving the cross sections 
constant. Nevertheless, also the cross sections might change in the medium. For the process 
NN— »NAK, the NAK vertex has to be considered which vanishes on the mean field level. 
Hence, changes of the cross section are here of higher order. They have to be computed 
by taking into account the p-wave interactions of nucleons and kaons and will change the 
angular distribution of the produced kaons in heavy ion reactions. The investigation of these 
effects is an interesting task and will be considered in a forthcoming work. 
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TABLES 

TABLE I. The coupling constants of the parameter sets used. The vector coupling constant for 
the A are taken from SU(6)-relations. The coupling constants of the kaons to the a- and 5-meson 
are fixed by the s-wave KN-scattering lengths. The vector coupling constants are chosen from 
SU(3)-relations. The parameters for the scalar and vector selfinteraction terms are not given, they 
can be found in the corresponding references. 
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FIGURES 
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FIG. 1. The Schrodinger equivalent potential of the A for several parameter sets of the RMF 
model as a function of density. The curve labelled MDG is the non-relativistic potential fit to 
hypernuclear data of Dover, Millener and Gal pif| . 
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FIG. 2. The energy of kaons and antikaons in nuclear matter as function of density for the 
soft EOS (parameter set TM1). 
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FIG. 3. The same as Fig. |2| for the hard EOS (parameter set NL-Z). 
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FIG. 4. The Q-values of the production processes of kaons and antikaons (the two upper 
bunches of curves) and the annihilation processes of antikaons and A's (lower bunches of curves) 
versus the density for the soft EOS. 
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FIG. 5. The same as Fig. | for the hard EOS. 
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